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ABSTRACT 

Aims. We investigate the dust-to-gas mass ratio and the environmental effects on the various components of the interstellar medium 
for a spatially resolved sample of Virgo spirals. 

Methods. We have used the IRAM-30m telescope to map over their full extent NGC 4189, NGC 4298, NGC 4388, and NGC 4299 
in the I2 CO(1-0) and the !2 CO(2-l) lines. We observed the same lines in selected regions of NGC 4351, NGC 4294, and NGC 4424. 
The CO observations are combined with Herschel maps in 5 bands between 100-500 from the HeViCS survey, and with HI data 
from the VIVA survey, to obtain spatially resolved dust and gas distributions. We studied the environmental dependencies by adding 
to our sample eight galaxies with I2 CO(1-0) maps from the literature. 

Results. We estimate the integrated mass of molecular hydrogen for the galaxies observed in the CO lines. We find molecular-to-total 
gas mass fractions between 0.04 < f mol < 0.65, with the lowest values for the dimmest galaxy in the B-band. The integrated dust-to-gas 
ratio ranges between 0.01 1 and 0.004. For the 12 mapped galaxies we derive the radial distributions of the atomic gas, molecular gas, 
and dust. We also study the effect of different CO-to-H 2 conversion factors. Both the molecular gas and the dust distributions show 
steeper radial profiles for Hi-deficient galaxies and the average dust-to-gas ratio for these galaxies increases or stays radially constant. 
On scales of ~ 3 kpc, we find a strong correlation between the molecular gas and the 250 jim surface brightness that is tighter than 
average for non-deficient galaxies. The correlation becomes linear if we consider the total gas surface mass density. However, the 
inclusion of atomic hydrogen does not improve the statistical significance of the correlation. 

Conclusions. The environment can modify the distributions of molecules and dust within a galaxy, although these components are 
more tightly bound than the atomic gas. 

Key words. Galaxies: evolution - Galaxies: clusters: individual: Virgo cluster - Galaxies: individual: NGC 4189, NGC 4294, NGC 
4298, NGC 4299, NGC 4351, NGC 4388, NGC 4424 - ISM: Molecules - Infrared: dust 



1. Introduction 

The stellar and gas distributions in a cluster galaxy can be drasti- 
cally modified through environmental effects. In Virgo, a young, 
close id ~ 17 Mpc) and still dynamically active cluster (see 
Bose lli & Gavazzi} |2006| and references therein), the environ- 
ment effects are clearly seen for the atomic gas. Through tidal 
interaction and/or ram pressure stripping a spiral galaxy can lose 

* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 

** Based on observations carried out with the IRAM 30m Telescope. 
IRAM is supported by INSU/CNRS (France), MPG (Germany) and 
IGN (Spain). 



a considerable fraction of its neutral gas and become deficient in 
HI (|Davies & Lewis||1973||Giovanardi et al. 1983; Haynes et 



[aT] |1984[|Cayatte et al.||1990||Chung et al.| j2009 ) 



For the molecular gas-phase the situation is more complex. 



The typical density of a molecular cloud core is 10 cm 



-3 



( Caselli et al. 2002), and the molecular gas is more tightly bound 
in the galaxy potential, being confined to the inner part of the 
disk which may prevent gas removal via stripping. According 
to |Kenney & Young] ( [T9 86 1989), Virgo spirals show no clear 
signs of deficiency in the molecular component. In the Coma 
cluster, where the ram pressure is expected to be stronger due 
to the high ICM densities, there is no clear evidence of molecu- 
lar deficiency (Casoli et al. 1991 . However, as noted by Boselli 
|& GavazzT| ( |2006| l, the sample of Casoli et al7| ( |1991[ ) is selected 
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according to IR luminosity, and their results are biased toward 
strong CO emitters (see also|Boselli et al.| 1995a 1997 1. Using 



the same data of Kenney & Young (1989 ), Rengarajan & Iyengar 
( 199 2| found evidenc e of a CO deficit in some Virgo spiral, and 



Voll mer et aL| ( |2008) l detected in NGC 4522 molecular gas dis 
placed from the galactic disk because of ram pressure. 



Boselli et al. (2002i, |Fumagalli & Gavazzi (2008 1 and 



Fumagalli et al. (2009 1 found a molecular depletion in highly HI 
deficient galaxies, where the HI is depleted well inside the opti- 
cal disk. Since not all the Hi-deficient galaxies examined in their 
papers are cluster galaxies, it is of interest to analyze a sample of 
normal and Hi-deficient galaxies in a cluster. Also, it is still not 
clear if the depletion of the molecular component in Hi-deficient 
galaxies is due to ram pressure stripping of the molecular com- 
ponent, or to the lack of HI from which molecular clouds form. 
T he Herschel Virgo Cluster Survey (HeViCS; |Davies et al.| 



2010 



Auld et al. 



2012 1 has mapped an area of 84 square deg z in 



the Virgo cluster, and one of the goals of the project is to investi- 



gate mechanisms of gas recycling in clusters of galaxies. Cortese 



et al. (20121 analyzed the galaxies in the Herschel Reference 
Survey (HRS; Boseili et al. 12010), including many Virgo galax- 



ies observed within HeViCS, and found that the global dust-to- 
stellar mass ratio decreases as a function of HI deficiency, pro- 
viding strong evidence of dust stripping. Furthermore they found 
that at fixed stellar mass, the dust-to-HI mass ratio decreases 
with the HI deficiency, with a correlation significantly weaker 
than for the dust-to-stellar mass ratio. They argue that this mild 
correlation is a consequence of the greater extension of the HI 
disks which produces a rapid decline of the dust-to-HI ratio in 
the outer regions of galaxies. For a sample of 35 HeViCS Virgo 
spirals with available integrated CO fluxes, Corbelli et al. ( 2012| > 
confirmed that the dust-to-stellar mass ratio decreases and the 
dust-to-total gas ratio increases when the galaxies are more HI- 
deficient. 

The investigation of the internal properties of galaxies can 
help understand environmental effects. Cortese et al. (20101 in- 
vestigated 15 Virgo spiral galaxies, and found that the radial dis- 
tributions of dust in Hi-deficient galaxies are truncated, already 
in the inner regions of the disks. 

It is therefore important to study the relation between the ra- 
dial distribution of dust and the different gas components within 
a spiral disk. To reach this goal, it is necessary to have a sample 
of galaxies observed in all the interstellar medium (ISM) compo- 
nents. The proximity of Virgo allows a good resolution of the gas 
and dust components (a typical beam size of 20" corresponds 
to a resolution of 1.6 kpc). Several Virgo spiral galaxies have 
been mapped with the Very Large Array (VLA) in HI, within the 
survey "VLA Imaging of Virgo spirals in Atomic gas" (VIVA; 
Chung et al.||2009"| . About half of them are mapped in the FIR 
within the HeViCS project ( |Davies et al.| |2010[ ). However, only 
a subset of these have observations that trace the spatial distri- 
bution of the molecular component, major axis profiles or maps 
such as those shown in |Chung et al.| ( [2007| l, |Kuno et al.| ( |2007| l. 
To complete the sample we thus began a project to map the CO 
emission in the remaining low luminosity objects. 

In this paper we present 12 CO(1-0) and 12 CO(2-l) emission- 
line observations for seven galaxies, obtained at the IRAM-30m 
telescope. We use these new data to investigate the balance be- 
tween different gas phases in spiral disks and their relation with 
the dust component as a function of the environment. The paper 
is organized as follows: in Sec. [2] we describe the sample, and 
the ancillary data used in the paper. In Sec.[3]we present the new 
IRAM-30m observations and the relative results. In Sec. [4] we 
discuss in detail the radial distribution of the dust and gas phases 



for the galaxies mapped in our new campaign and for additional 
Virgo spirals with available CO, HI, and dust maps. We also in- 
vestigate the variations of the dust-to-gas ratio across the disk. 
In Sec. [5] we perform a pixel-by-pixel analysis of the relation be- 
tween the far-infrared (FIR) emission and the gas surface mass 
density. The main conclusions are summarized in Sec. [6] 



2. The sample and ancillary data 

We have compiled a sample of bright (mg < 13) late-type galax- 
ies in the Virgo cluster with available VIVA and HeViCS ob- 
servations. Here we present 12 CO(1-0) and 12 CO(2-l) emission- 
line maps of four spirals (NGC 4189, NGC 4298, NGC 4299, 
NGC 4388) obtained with the IRAM-30m telescope. We also 
observed a few selected regions of 3 additional low luminos- 
ity galaxies, NGC 4351, NGC 4294, and NGC 4424. To enlarge 
the sample of fully mapped galaxies we also consider eight ad- 
ditional galaxies mapped in 12 CO(1-0) at the Nobeyama 45 m 
telescope by |Kuno et al. ( 2007) at comparable spatial resolu- 
tion. This sample of 12 mapped galaxies spans a wide range of 
morphological types, from Sa (NGC 4424) to Scd (NGC 4299). 
Since we are interested in studying the environmental effects on 
the gas distribution, we excluded from the sample gravitationally 
interacting pairs. 

Interferometric observations of HI for the VIVA sample are 
described in Chung et al. (2009 1. The velocity resolution of 
VIVA (10 km s" 1 ) is comparable to the final resolution of our 
CO data (10.4 km s~'). The sensitivity of the VIVA observa- 
tions is around 3-5 x 10 19 HI atoms cm -2 (3cr, 0.3 M Q pc" 2 ) 
at a typical angular resolution of ~25". We quantify the HI de- 
ficiency, defm, m our sample using the prescription of Chung 
|et al.| ( |2009] ), which assumes a mean HI mass-diameter relation, 
independently of morphological type. Our galaxies span a wide 
range of values, -0.43 < defm ^ 1.47, revealing different stages 
of interaction with the environment. Even if we adopt an HI de- 



ficiency criterion based on morphological type (e.g. |Giovanardi 
et al. 1983| l, our sample still covers a wide range of defm- from 
-0.2 to 1.1. Given that the uncertainties in defm are ~ 0.25-0.3 
dex we shall consider a galaxy as non-deficient only if it has 
defm < 0.4. 

FIR images for dust emission at 100, 160, 250, 350, and 500 



fim come from the HeViCS full depth data set (see Auld et al. 
2012, for details on the data reduction). The beam Full Width 
Half Maximum for these wavelengths are 9", 13", 18'.'2, 24"5, 
and 36", respectively. 



The dust surface mass density is derived as in Magrini et al 
20T"fl i (see also |Smith et al.| [20T0l |Bendo et al.| |2010a| >. The 



five images in the Herschel photometric bands were convolved 
to the poorer resolution of the 500 /vm, then registered on the 
same pixel grid. Each pixel was fitted with a single temperature 
modified black body, with dust emissivity k oc \P and /3 = 2; 
the filter response function was taken into account for the color 



correction. As shown by Davies et al. ( 201 l| l for global fluxes, 
and by Magrini et alT]( |201 1) for pixel-by-pixel surface brightness 
analysis, our choice of the dust emissivity provides a satisfactory 
fit to the spectral energy distribution. If we exclude the notorious 
uncertainty of the dust emissivity, the uncertainties due to the 
fitting procedure amounts to 10-15% for the full depth HeViCS 
data set. Other sources of uncertainty in the determination of the 
total mass are discussed in Magrini et aL]( |201 l| l. 

In Table [T] we summarize the main properties of the galaxies 
in the sample, and in Appendix|A]we describe more in detail the 
galaxies for which we present new CO observations. 
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Table 1. General properties of the sample 



Galaxy 


RA 


Dec 


Type 


m B 


D 


Vkel 


a 


i 


P.A. 


def m 


Z, 




deg 


deg 






Mpc 


km s" 1 




deg 


deg 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 








Targets of our observations 












NGC4189 


183.447 


13.42 


Sc 


12.69 


32 


2114 


2.40 


45 


66 


0.25 


9.09 


NGC 4294 


185.324 


11.51 


Sc 


12.68 


17 


355 


3.24 


70 


155 


-0.11 


8.95 


NGC 4298 


185.386 


14.61 


Sc 


11.95 


17 


1136 


3.24 


57 


140 


0.41 


9.08 


NGC 4299 


185.419 


11.50 


Scd 


12.99 


17 


232 


1.74 


22 


42 


-0.43 


8.83 


NGC 4351 


186.006 


12.20 


Sc 


12.99 


17 


2316 


1.20 


49 


61 


0.23 


8.88 


NGC 4388 


186.445 


12.66 


Sab 


11.87 


17 


2515 


5.62 


75" 


92 


1.16 


9.09 


NGC 4424 


186.798 


9.42 


Sa 


12.33 


23 


437 


3.63 


62 


95 


0.97 


9.08 



Galaxies with 12 CO(1-0) maps from Kuno et al. (20071 



NGC 4192 


183.451 


14.90 


Sb 


10.73 


17 


-139 


9.77 


78 


155 


0.51 


9.08 


NGC 4254 


184.707 


14.42 


Sc 


10.45 


17 


2404 


5.37 


30 


68 


-0.10 


9.09 


NGC 4321 


185.728 


15.82 


Sc 


10.02 


17 


1575 


7.41 


33 


30 


0.35 


9.07 


NGC 4402 


186.532 


13.11 


Sc 


12.64 


17 


230 


3.89 


78 


90 


0.74 


9.08 


NGC 4501 


187.996 


14.42 


Sbc 


10.50 


17 


2282 


6.92 


59 


140 


0.58 


9.05 


NGC 4535 


188.585 


8.20 


Sc 


10.73 


17 


1964 


7.08 


46 





0.41 


9.09 


NGC 4569 


189.208 


13.16 


Sab 


10.08 


17 


-221 


9.55 


65 


23 


1.47 


9.08 


NGC 4579 


189.431 


11.82 


Sab 


10.52 


17 


1516 


5.89 


38 


95 


0.95 


9.06 



Column 1: name; col. 2-3: equatorial coordinate (J2000); col. 4: morphological type (Binggeli et al. 1985 1993); col. 5: B-band magnitude; col. 
6: distance; col. 7: heliocentric velocity; col. 8: major axis optical diameter; col. 9: inclination; col. 10: Position Angle; col. 11: HI deficiency 
parameter; col. 12: central oxygen abundance expressed as 12+log O/H from the mass-metallicity relation of |Trem onti et al. ( 2004 1. Data sources: 
col. 1-7, GOLDMine l |Gavazzi et aLl [2003 ); col. 8-11, [Chung etal] p00"9} ; col. 12, |Corbelli eTaIJ l [20l2l >. 
8 The inclination of this galaxy has been computed directly from the axis ratio. 



3. IRAM 30-m observations 

The CO observations were carried out at the IRAM 30-meter 
millimeter-wave telescope on Pico Veleta (Spain) during two ob- 
serving runs in June 2010 and April 201 1. 



3.1. Observational setup and data reduction 

We observed the 12 CO(1-0) and 12 CO(2-l) rotational transitions, 
with rest frequencies of 1 15.271 and 230.538 GHz, respectively. 
As front-end we used the Eight Mixer Receiver (EMIR; |Carter| 
et al. 2012 1, and adopted three different back-ends: 



- the 4MHz filterbank, used for the 12 CO(1-0) line only, pro- 
viding a total band width of 4GHz, and a channel width of 
4MHz, i.e. 10.4 km s" 1 ; 

- the Wideband Line Multiple Autocorrelator (WILMA) with 
a bandwidth of 1860 channels (4GHz), and a spectral reso- 
lution of 2MHz. This corresponds to a spectral resolution of 
5.2, and 2.25 km s 1 for the 12 CO(1-0) and 12 CO(2-l) lines; 

- the VErsatile SPectrometer Assembly (VESPA) for which 
we selected the resolution of 1 .25MHz, with a smaller band- 
width, 420 MHz. This setup gives a spectral resolution of 
3.25, and 1.62 km s 1 for the 12 CO(1-0) and 12 CO(2-l) lines; 

We performed a flux calibration every 15 minutes, checked 
the pointing every 2h, and the focus every 6h, as customary. The 
spectra were calibrated in main-beam temperature, assuming the 
main-beam and forward efficiencies 77,,,/, = 0.78 and rjf or = 0.95 
for the 12 CO(1-0) and j] mb = 0.59, r] for = 0.91 for 12 CO(2-l), 
according to the values given in the IRAM-30m website as of 
July 2011. The main-beam temperature is T m t, = T ]f or /r] m b- T* A , 
where T* is the antenna temperature ( |Wilson et al. |2009] >. 

In order to increase the S/N ratio we averaged the spectra 
taken with all the backends to a common resolution of 10.4 km 



s 1 . We removed the baseline outside the signal window by fit- 
ting polynomial functions^] of orders between 1 and 3. The half 
power beam widths (HPBW) are 22/5 for 12 CO(1-0) and 11725 
for 12 CO(2-l). We observed our sample in 2 different observing 
modes: On The Fly (OTF) map and Position Switching (PS). The 
choice of the observing mode was based on the estimated total 
CO brightness from FIR emission. 

In OTF mode the telescope moves smoothly across the tar- 
get and collects spectra while scanning. We set the integration 
time per position to 2 sec and the distance between each posi- 
tion along the scan to 4", below the critical Nyquist frequency. 
After each scan a reference spectrum is taken in a region with 
no sources. Several scans along a chosen direction were taken to 
cover an area of about 1 square arcmin. The observations were 
repeated several time along parallel and perpendicular scan di- 
rections to reach the required sensitivity. In total, we observed 
NGC 4189 for 8.1h, NGC 4298 for 4h, NGC 4388 for 2.4h, and 
NGC 4299 for 4h. The sensitivities we achieved for the 12 CO(l- 
0) and 12 CO(2-l) lines are shown in columns 3, 4 of Table [2] 
The size of each map is given in column 2 of Table [2] Scans 
were done along the RA and DEC directions, with the exception 
of NGC 4298, where a scan direction parallel and perpendicular 
to the major axis was used. 

The map pixel size is chosen in order to achieve roughly 
Nyquist sampling, 973 for 12 CO(1-0) and 477 for 12 CO(2-l) 
(HPBW/2.4). Spectra at each pixel are summed using the sigma- 
weighted average to obtain the map. The observations have been 
smoothed with a gaussian of size 1/3 of the telescope beam to 
gain i n sensitivity at the expense of a slightly larger beam size 
(see|Mangum et al. 2007 1. In the following analysis, the 12 CO(2- 



1 Most of data reduction was performed using GILDAS software: 
http://www. imm.fr/IRAMFR/GlLDAS 
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Table 2. Integrated quantities for observations in OTF mode 



















M l-0 


Galaxy 


Map Size 


cr i-o 




Ji-0 


h-i 


h-U'i-o 


y H 2 




" x " 


mK 


mK 


Kkms- 1 


Kkms 1 




M pc- 2 


10 s M Q 


(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


NGC4189 


96 x 75 


7 


10 


3.73±0.34 


3.43±0.78 


0.94 


11.95±1.1 


17.4+1.3 


NGC 4298 


120 x 60 


11 


14 


19±0.66 


6.04±2.4 


0.76 


23.32±2.1 


10.4±0.8 


NGC 4299* 


104x104 


19 


42 


1.28±0.28 


1.2±0.7 




4.1±2.6 


1.1+0.2 


NGC 4388 


117x 42 


13 


39 


8.93±0.75 


10.8±2.6 


0.9 


28.55±2.4 


7.13±0.6 



Column 1: name; col. 2: map size; col. 3, 4: measured uncertainties in 10.4 km s _I channels for a 22'.'5 beam; col. 5, 6: average value of the 
integrated intensity; col. 7: average 12 CO(2-1)/ 12 CO(1-0) intensity ratio; col. 8: average surface mass density of H 2 using 12 CO(1-0) and the 
Galactic Xqo', col. 9: total IF mass. 

*For this galaxy, because of the low S/N ratio, the results shown were obtained integrating inside an aperture of R = 35" (see Sec. |3.2| for details). 
Table 3. Integrated quantities for observations in PS mode 



Galaxy 


ARA, ADec 


Vl-0 


02-1 


Ii-o 


I2-1 


h-i/h-o 


yH 2 


M Hl 






mK 


mK 


Kkms 4 


K km s 1 




M Q pc - 


10 8 M 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


NGC 4351 


+0.0,+0.0 


4 


13 


2.29+0.14 


2.3+0.31 


1.00 


7.34+0.44 


0.71 


NGC 4294 


+0.0,+0.0 


4 


4 


2.66+0.11 


1.36+0.12 


0.51 


8.5+0.35 


0.73 




+9.3,-19.9 


4 


8 


1.08+0.18 


0.99+0.23 


0.92 


3.47+0.58 




NGC 4424 


+0.0,+0.0 


6 


10 


9.41+0.24 


9.0+0.4 


0.96 


30.1+0.77 


4.2 




+21.9,-1.9 


3 


13 


2.28+0.16 


2.7+0.5 


1.18 


7.31+0.52 






+43.8,-3.8 


11 


11 


1.08+0.25 


1.02+0.37 


0.94 


3.67+0.8 






-21.9,+ 1.9 


11 


13 


6.31+0.4 


5.7+0.49 


0.9 


20.2+1.28 






-43.8,+3.8 


11 


11 


1.05+0.16 


1.04+0.5 


0.99 


3.35+0.52 





Column 1: galaxy name; col. 2: observed offset with respect to the galaxy center; col. 3, 4: lcr rms associated to 12 CO(1-0) and 12 CO(2-l) lines, 
respectively; col. 5, 6: integrated CO flux (/ Tdv) for 12 CO(1-0) and 12 CO(2-l); col. 7: 12 CO(2-l)/ 12 CO(1-0) intensity ratio; col. 8: surface mass 
density derived from the 12 CO(1-0) line using the Galactic Xco; col. 10: total H2 mass (see Sec. 3.2 for details). 



1 ) maps have been convolved with an appropriate bidimensional 
Gaussian kernel to match the resolution of the 12 CO(1-0). 

To estimate the Moment and Moment 1 maps we defined in 
each spatial pixel a spectral signal window and integrated within 
it. We calculated the noise <x outside the window and then esti- 
mated the noise related to the signal as N = cr ■ Av / ~\/N c i,, where 
Av is the width of the window in km s _1 and N c h is the number 
of channels inside the signal window. In the final cube we took 
into account only the pixels with S/N > 3. 

Three galaxies with expected fainter CO emission were ob- 
served in PS mode at a few selected positions, with the telescope 
switching between the galaxy and an offset position for back- 
ground subtraction. In these observations we integrated for 15 
minutes on each position and achieved the sensitivities shown in 
col. 3 and 4 in Table [3] The original plan was to observe several 
positions along the major axis spaced by the I2 CO(1-0) beam- 
size. However, due to bad weather conditions and time limita- 
tions this was possible only for NGC 4424, which was observed 
in 5 positions. For NGC 4294 it was possible to observe only the 
galaxy center and another position, and finally for NGC 435 1 we 
only observed the galaxy center. The observed positions with re- 
spect to the center are shown in col. 2 of Table [3] 

3.2. Results on individual galaxies 

Table [2] shows the average value of the integrated intensity for 
12 CO(1-0) and 12 CO(2-l) lines (column 5 and 6), together with 



the average intensity line ratio (column 7) for galaxies observed 
in OTF mode. In column 8 and 9 we report the average sur- 
face brightness and the total mass of molecular hydrogen, ob- 
tained integrating over the pixels with surface mass densities 
above 3cr. We converted the 12 CO(1-0) integrated intensities 
into molecular gas masses using the standard Galactic conver- 
sion factor X co = N H JI C o = 2 • 10 20 cirT 2 [K km s -1 ] -1 , esti- 



mated in the solar neighborhood ( Strong & Mattox 1996 1 Dame 
|et al.||200T] l. In Sec.|4]we investigate how the choice of different 
Xqo factors affects the estimate of molecular hydrogen surface 
mass densities. We estimate the molecular fraction, defined as 
fmoi - M Hi /(Mh 2 + M HI ), using the mass of atomic hydrogen 
given by Chung et al. ( 2009| l. 

In Figs.[T]|3]|4]the CO maps obtained in OTF mode are pre- 
sented in the following order: in the top left and right panels we 
show the moment and moment 1 map of the 12 CO(1-0) inten- 
sity emission; in the middle left and right panels the moment 
12 CO(2-l) map and the 12 CO(2-1)/ 12 CO(1-0) line intensity ratio 
are shown. The channel velocity maps with Av = 10.4 km s _1 
are given in the bottom panels. 

NGC 4189 

The CO maps of NGC 4189 are shown in Fig. [T] The 3cr noise 
in 12 CO(1-0) and 12 CO(2-l) maps is 4.2 and 4 M G pc~ 2 , respec- 
tively. The total mass of molecular gas is Mh 2 = 17.4 ± 1.3 • 10 8 
M Q , and the molecular-to-total gas mass fraction is /„,„/ = 0.43. 
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Fig. 1. NGC 4189 maps. Top left: 12 CO(1-0) Moment map, contours: 4.2, 6.9, 9.6, 12.3, 15, 17.7 M Q pc -2 . Top right: 12 CO(1-0) 
Moment 1 with 2000 < v < 2230 and Av = 28 km s _1 . Middle left: 12 CO(2-l) Moment map, contours: 4, 7.25, 10.5, 13.75, 17, 
20.25, 23.5 M G pc~ 2 . Middle right: 12 CO(2-1)/ I2 CO(1-0) ratio. The beam size is illustrated in the bottom left of each panel. Bottom 
panel: channels map of 12 CO(1-0) with 1980 < v < 2250 km s" 1 and Av = 10 km s" 1 . 
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Fig. 2. Ha emission (from GOLDMine) of NGC 4189 (top row) and NGC 4298 (bottom row) with overplotted 12 CO(1-0) (left 
column), 250 /mi emission (middle column), and HI surface mass density (right column). The contours of 12 CO(1-0) line are 4.2, 
6.9, 9.6, 12.3, 15, and 17.7 M G pc~ 2 for NGC 4189 and 0.5, 8.8, 17.1, 25.4, 33.7, 42, and 50.3 M pc 2 for NGC 4298. The 250 /mi 
contours are at 0.23, 0.28, 0.31, 0.36, 0.42, 0.48, 0.56, and 0.65 Jy beanT 1 for NGC 4189 and 0.2, 0.3, 0.41, 0.63, 0.84, and 1.06 Jy 
beam 1 for NGC 4298. The HI surface mass density contours are at 4.8, 7.2, 10, 12.4, 14.8, and 17.2 M G pc~ 2 for NGC 4189, and 
at 4.3, 5.2, 6, 7, 8, 9.1, 10.6 M Q pc~ 2 for NGC 4298. The beam size of each observations is shown in the bottom left of each panel. 



The 12 CO(1-0) emission is less extended (~ V.6 at 4.2 M Q pc~ 2 ) 
than the atomic hydrogen (~ 2'.6 at 4.8 M pc~ 2 ) but comparable 
to the extent of the optical disk (2!26). The total intensity of the 
12 CO(2-l) lines is comparable to that of the 12 CO(1-0) line. 

The 12 CO(1-0) emission (top left panel of Fig. [TJ has a reg- 
ular structure centered in the South-East region and is more ho- 
mogeneously distributed with respect to the I2 CO(2-l) emission, 
which has a clumpier structure dominated by a few unresolved 
regions (of 1 pixel size). Taking into account only regions with 
S/N > 3 the 12 CO(2-1)/ 12 CO(1-0) ratio varies between 0.4 and 1 
(bottom right panel of Fig. [TJ. These values are consistent with 
the bulk of the emission being optically thick. 

In the top row of Fig. [2] we compare qualitatively the Ha 
emission of NGC 4189 (from GOLDMine; |Gavazzi et aT][2003] l 
with the 12 CO(1-0) (left), HI (middle), and the 250 fan (right) 
emission. Dust emission at 250 /mi is comparable in size with the 
12 CO(1-0) emission. The distribution of dust emission across the 
disk is regular and more similar to the CO brightness distribution 
than to the HI. There is a higher concentration of the dust in the 
South-East regions, in agreement with a similar feature seen in 
the HI, 12 CO(1-0), and Ha maps. 

NGC 4298 

The CO maps for NGC 4298 are shown in Fig. [3] We estimate 
a total mass M H , = 10.4+0.8 • 10 s M and f mo , = 0.65. Both 



the 12 CO(1-0) (top left panel of Fig.[3j» and 12 CO(2-l) emissions 
(middle left panel of Fig. [3} peak in the central regions, con- 
sistently with the observed Ha emission (bottom row of Fig. [2] 
from GOLDMine). 

The Moment 1 map (top right panel of Fig. [3]l shows a quite 
regular structure in the northern part of the disk and slightly dis- 
turbed velocity contours in the South. The 12 CO(2-l)/ 12 CO(l- 
0) ratio has values between 0.3 and 1.4, with the highest val- 
ues found in isolated clumps of 1-2 pixels in size. The location 
of these pixels corresponds to the brightest HII regions in the 
galaxy. 

Both the HI and the CO emission have a higher surface mass 
density in the South-East region, but the 250 /mi emission is only 
slightly asymmetric (see bottom row panels of Fig. |2]». 

NGC 4388 

For NGC 4388 the CO maps are shown in Fig. [4] The estimated 
total mass of molecular hydrogen is 7.1 ±0.6 ■ 10 8 M Q and the 
molecular fraction is f mo \ = 0.63. Unfortunately the I2 CO(2-l) 
emission has SNR > 3 only in the central regions and the relative 
map is smaller than 12 CO(1-0) map. Both the HI and CO distri- 
butions throughout the disk drop sharply below the sensitivity of 
our observations; this can be seen in the middle panel of Fig. [5] 
where we overlay the contours of the molecular gas onto a refer- 
ence g-band image from SDSS-DR6 ( |Adelman-McCarthy et al.| 
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Fig.3. NGC 4298 maps. Top left: 12 CO(1-0) Moment map, contours: 0.5, 8.8, 17.1, 25.4, 33.7, 42, and 50.3 M Q pc~ 2 . Top right: 
12 CO(1-0) Moment 1 with 1035 < v < 1235 and Av = 20 km s _1 . Middle left: 12 CO(2-l) Moment map, contours: 0.5, 7.6, 14.8, 
21.9, 29, and 36.2 Mq pc~ 2 . Middle right: 12 CO(2-l) over 12 CO(1-0) ratio. The beam size is shown in the bottom left of each panel. 
Bottom panel: channels map of 12 CO(1-0) with 985 < v < 1275 km s _1 and Av = 10 km s _1 . 
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Fig. 4. NGC 4388 maps. Top left panel: 12 CO(1-0) Moment map, contours: 1, 12, 23, 34, 45, 56, and 67 M Q pc -2 . Top right panel: 
12 CO(1-0) Moment 1 with 2370 < v < 2800 and Av = 43 km s _1 . Middle left panel: 12 CO(2-l) Moment map, contours: 0.01, 8.6, 
17.3, 26, 34.7, 43.3, and 52 M Q pc~ 2 . Middle right panel: 12 CO(2-l) over 12 CO(1-0) ratio. The beam size is shown in the bottom 
left of each panel. Bottom panel: channels map of 12 CO(1-0) with 2304 < v < 2744 km s _1 and Av = 10 km s _1 . 
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Fig. 5. Top panel: 12 CO(1-0) em ission of NGC 4388 with over- 
plotted HI emission from VIVA ( |Chung et al.| [2009 ). Contours: 
15.7, 18.7, 22.1, 25.2, 28.6 M Q pc z . Middle panel: SDSS g-band 
image of NGC 4388 with overplotted 12 CO(1-0) emission, con- 
tours: 1, 12, 23, 34, 45, 56, and 67 M pc~ 2 . Bottom panel: SDSS 
g-band image of NGC 4388 with overplotted 250 fim dust emis- 
sion, contours: 0.35, 0.71, 1.07, 1.43, 1.79, and 2.15 Jy beam- 1 . 
The beam size is shown in each panel. The arrow in the bottom 
panel indicates the direction of ram pressure. 

NGC 4299 Spectrum R = 35" 




Fig. 6. Letf panel: Ha emission of NGC 4299 with overplotted a 
circle with R = 35" selected to integrate the 12 CO(1-0) emission. 
Right panel: averaged spectra of NGC 4299 inside the circle de- 
fined in left panel. The vertical shaded lines show the velocity 
limits used for the integration. 



2008 ). The atomic and molecular components are truncated at 



R/Ris — 0.6, where R25 is the radius at a B-band surface bright- 
ness of 25 mag arcsec" 2 (R25 — a/2, see column 8 of Table fTb, 
as expected for highly Hi-deficient galaxies {def m = 1.16). The 
top panel of Fig.[5]shows that the radial extent of the atomic and 
molecular phases for this galaxy are comparable. 

The structure of the 12 CO(1-0) emission is asymmetric, with 
a more prominent elongated diffuse tail in the east side of the 
disk (middle panel of Fig. |5). A similar feature is present also 



in the HI distribution ( |Chung et al. 2009[ ), and this overden- 



sity is higher where the gas plume observed in |Oosterloo & van 
|Gorkom 



(2005 ) connects with the disk. At 250 jum the disk is 
also slightly asymmetric, similar in shape to the CO distribution 
(middle and bottom panels of Fig . [5]l . This asymmetry in the ISM 
component is probably due to a recent ram pressure episode (see 





v (km/s) 



Fig. 7. Top panel: SDSS g-band image of NGC 4351. White cir- 
cles show the beam size for 12 CO(1-0) (solid line) and 12 CO(2- 
1) (dashed line). Bottom panel: 12 CO(1-0) (black solid line) and 
I2 CO(2-l) (red dashed line) observations in the position showed 
in the left panel. 



arrow in Fig.[5] [Vbllmer & Hu chtmeier 2003; Pappalardo et al. 
2010). Because of the spatial limitation of our map, we were un- 
able to investigate whether there is any stripped molecular gas 
associated with the elongated HI tail seen by Oosterloo & van 
Gorkom| (2005 1. The channel maps (top right panel of Fig. [4]) 
do not show the typical "spider diagram" of spiral galaxies, be- 
cause NGC 4388 is highly inclined (i = 75°) and the isovelocity 
contours are crowded along the line of sight. 

NGC 4299 

Because of the low signal-to-noise ratio of the observations, the 
CO emission in this galaxy has been determined by integrating 
the spectra in the velocity range 130 < v < 300 km s _1 as deter- 
mined from the HI emission ( |Chung et al.| |2009| ), and inside a 
circle of 35" of radius from the galaxy center (Fig.|6|. The size 
of the circle was chosen to match the Ha emission, taking into 
account all the star forming regions of the disk. The integration 
gives an intensity Iqo = 1-28 K km s _1 , with a noise of cr = 0.28 
K km s -1 (i.e. SNR ~ 4.5). The estimated mass inside the circle 
is Mh 2 — 1.1 ± 0.2 • 10 s Mq with an average surface mass den- 
sity of 4 M Q pc~ 2 and one of the lowest molecular mass fractions 
f mol = 0.08 (Tableg). 

NGC 4294, NGC 4351 and NGC 4424 

For the galaxies observed in PS mode Table [3] reports the inte- 
grated CO line emission (col. 5, 6), the line ratio (col. 7), and the 
surface mass density (col. 8) for each position. The total M#, 
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mass (col. 9) was estimated assuming an exponentially decreas- 
ing CO-intensity for the gas (Kenney et al. 1996 ). For both NGC 
4294 and NGC 4424 we estimated a scalelength of 0.1 R 2 $. Since 
NGC 435 1 was observed only in the central region, we assumed 
a scalelength 0.2 R25, as determined in Schruba et al. (201 1 1 for 
a sample of 33 nearby spiral galaxies. 

Fig. [7] [8] and [9] show the spectra with the associated posi- 
tions. For NGC 4351 and NGC 4294 we found the lowest H 2 
masses of the observed sample, 7.1T0 7 M G , and 7.3T0 7 M , and 
also low molecular fractions, f mo i = 0.17, and 0.04, respectively. 
For NGC 4424 we found a M Hl = 4.2 TO 8 M G , with a molecular 
fraction f mo \ = 0.51. The integrated fluxes of both 12 CO(1-0) and 
12 CO(2-l) have comparable values. 

3.3. Comparison with previous molecular mass estimates 
Figure [10] shows the comparison of our measurement of the 



molecular gas mass with previous estimates (|Stark et al. 



Young et al. 1995 1 Boselli et al. 1995b; Smith & Madden 



1986 



1997 



Chung et al. 2009). We have rescaled the literature values to 
match our choice of Xqo and distance. Overall, there is a quite 
good agreement with previous observations. The only excep- 
tion is NGC 4388 for which |Stark et al.| ( [1986] ) reported M Hl 
= 2.2-1 Q 9 M P 



i , a value ~30% higher than our estimate. The beam 
size of their telescope is larger, 100", and they observed the disk 
in 5 different positions, with an offset of 1 '.5 in East-West and 
North-South directions. Surprisingly they found in the East and 
West side of the disk an emission comparable to the value in 
the center of the disk. Our map does not extend farther than 80" 
along the disk; however the emission decreases steeply already 
at ~30-40" from the center. Our molecular mass estimate agrees 
with the value of |Kenney & Y oung ( 1986) which have sampled 
similar regions of this galaxy as |Stark et ah] ( |1986[ ), but with a 
smaller beam. 

|Cortes et al.| (|2006) carried out interferometric observations 
of the lz CO(l-0)emission line of the central regions of NGC 
4424, and found a molecular gas mass of 1.6- 10 8 Mq. Given the 
larger area covered by our observations, this is compatible with 
the mass value we find in the central beam and with our higher 
total mass value. 

Our observations have a higher sensitivity and better spatial 
resolution than previous ones. For most of the galaxies, ct^q is 
between 4 and 1 1 mK. Moreover all the galaxies observed in 
OTF mapping mode have not been mapped previously. 



4. Environmental effects on the radial distribution 
of gas and dust 

In this section we study the influence of the environment on the 
surface mass density distributions for the molecular, atomic, and 
dust components. We consider the new mapped galaxies together 
with 8 additional galaxies observed at the Nobeyama 45-m tele- 



scope ( Kuno et al. 2007 ). The beam size of this set of obser- 



vations tFT5 7T Tcorriparable to the IRAM-30m telescope beam 
(2275 for I2 CO(1-0) line), but the typical rms noise level is 40- 
100 mK, a factor of ten higher than the level achieved by our 
observations. 

We divide our galaxies in three groups: non-deficient galax- 
ies with def m < 0.4 (NGC 4299, NGC 4254, NGC 4189, NGC 
4321), galaxies slightly Hi-deficient, with 0.4 < def m < 0.7 
(NGC 4298, NGC 4535, NGC 4192, NGC 4501), and HI- 
deficient galaxies, with def m > 0.7 (NGC 4402, NGC 4579, NGC 




i.O 8.5 9.0 

Log(MH2/M ) - literature 

Fig. 10. Comparison of our results with previous observations. 



References: Young et al. ( 1995[) ( blue triangles), |Chung et al. 
d2009)) (red sta r), |Stark et al.| ( |1986| ) (b lack circles), |Smith & 
Madden] \\991) (magenta squares), and |Boselli et al.| l |1995b) 
(green diamonds). Arrows represents upper or lower limits. 



4388, NGC 4569). We use this classification to characterize the 
environmental effects in our sample. 

We consider the total gas surface density as the sum of the 
atomic and molecular gas component, since we did not take into 
account the helium contribution. To scale the quoted surface 
mass density for helium, one should multiply the recovered val- 
ues by a factor ~ 1.36. We also neglected the ionized gas, which 
is expected to be a small fraction of the total gas mass. 

4. 1 . Radial profiles 

The radial profiles are obtained averaging surface mass density 
in elliptical annuli of width 0.15 R/R25, using the inclination and 
the major axis position angle given in columns 10, 1 1 of Table[T] 
The average value in each annulus is plotted when the total num- 
ber of pixels with surface mass density above the threshold value 
is at least 20% of the total number of pixels inside the annulus. 
Finally the radial profiles are corrected to face-on values. 

Fig. 1 1 shows the radial profiles of the gas and dust com- 



ponents. The panels are ordered according to the defm param- 
eter of each galaxy: the top, middle, and bottom rows refer to 
the non-deficient, slightly Hi-deficient, and Hi-deficient galax- 
ies, respectively. The open blue circles refer to the molecular 
surface mass density for a constant Galactic Xco value. The ra- 
dial profile of the atomic component is plotted using red star 
symbols. Since the beam size of the HI observations is similar to 
that of the CO observations, we kept all images at their original 
spatial resolution (~ 20-25"). The total gas profile is shown with 
black diamonds. The triangles connected by the green solid line 
show the dust surface mass density at the resolution of 36". 

As a general trend, the molecular gas has a higher surface 
mass density than the HI gas in the inner regions. In the external 
regions the mass surface density of H2 decreases, and the rela- 
tion between the total gas and the HI distribution depends on the 
environmental effects undergone by the galaxy. The dominance 
of the molecular gas in central regions is consistent, for some 
galaxies in our sample, with the results obtained in a sample of 
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Fig. 8. Same as Fig.^for the two positions observed in NGC4294. The RA and Dec offsets with respect to the center position are 
indicated in the top right of each spectrum panel. 



33 nearby galaxies by Bigiel et al. (2008). By analyzing the sur- 
face mass density of the different gas phases they recovered a 
threshold at which most of th e gas is converted into mole cular 
form, i.e. 9 M pc~ 2 (see also |Krumholz et al.| |2008| |2009[ for 
theoretical modelling of this effect). In Fig.|l l|the horizontal dot- 
ted line marks this value of surface mass density. In the central 
regions of NGC4535, NGC4501, and NGC4579 the molecu- 
lar phase is already dominant where the HI surface mass den- 
sity is only 1 M pc~ 2 . However, these are Hi-deficient galaxies 
and variations of the threshold value may occur. NGC 4299, the 
galaxy with the lowest metallicity in our sample, shows a lack 
of CO. In the central regions the HI dominates the total gas sur- 
face mass density, even if for this galaxy it is above the threshold 
value. 

The estimate of the molecular gas profile strongly depends 
on the assumption for the Xco factor. We have used so far the 
Galactic solar neighborhood value, but there is ample evidence 
that Xco varies in other galaxies as a function of physical condi 
tions (e.g. |Maloney & Blackj [T988] [Israeli 
2000|IBoselli et al.||2002HIsrael et al 1|2003 
et al.f |2008||Narayanan et al.||2005||Leroy et al.||201 l||Shetty et[ 
al.j |201 1| |Narayanan et al.j |2012[ |Bell et al.| |2007] k Thus, we 
consider an alternative conversion factor, which depends on the 
metallicity, and varies radially following the metallicity gradient 
along the disk: 



1997 



Israel 



Barone et al.j 
2005 1 Bolattol 



Mr 



10 20 / 



Wo, 



R 



R 



-Z c + 8.69, 



(1) 



25 



with Z c equals to the oxygen abundance in central regions, ex- 
pressed in the usual form 12+log O/H. Since we do not have Z c 
and gradients for all the galaxies in our sample we assume the 
central metallicity estimated using the mass metallicity relation 
(Tremonti et al. 2004) and reported in col. 12 of TablefT] We fix 
the solar value at Z = 8.69, as recently found by |Asplund et al.| 
(2009|, and assume a radial metallicity gradient oc -0.4 R/R25, 



as found by Magrini et al. ( 20 1 1 ] > for a sample of Virgo spirals. 
This slope should reproduce the effect of a metallicity gradient. 



The blue dashed lines in Fig. 1 1 represent the surface molec- 
ular mass density calculated according to Eq. [T] The choice of 
a radially varying Xco as a function of metallicity affects sig- 
nificantly the radial profiles and can reduce the molecular mass 
content up to 60%. Considering for our sample a typical metal- 
licity Z t . = 9.09 and according to Eq. [T| we obtain for the cen- 
tral regions Xco- 0-78 ■ 10 20 [K km s considerably lower 
than the Galactic X C o- For NGC 4189, if we use the galactic 
Xco we obtain a centrally peaked molecular distribution with a 



surface mass density 2-3 times larger than HI in the central re- 
gion. For R/R25 <; 0.7, when the surface mass density is below 
4 M Q pc~ 2 , the HI phase dominates over the molecular phase. 
Rescaling the conversion factor as in Eq. T[ because of the high 
value of Z c = 9.09, the atomic hydrogen becomes the dominant 
component already at R/R25 S 0.3, and the surface mass density 
of the molecular gas is never above 6 M Q pc~ 2 . For NGC 4298 
the molecular gas in the central region is always dominant for 
every choice of the conversion factor, and the effect of a rescal- 
ing is to reduce the surface mass density value according to the 
metallicity. For both NGC 4 189 and NGC 4298 the effect of a 
radially varying conversion factor is negligible in the external 
regions, R/R25 > 0.6, where, according to our assumed gradient, 
the metallicity approaches solar. For NGC 4388 the molecular 
component with a variable Xco has surface mass densities com- 
parable to the atomic component, as a consequence of the strong 
Hi-deficiency (def m = 1.16). For NGC 4299 (top left panel of 
Fig. 11 1 the effect on the total gas is negligible, because the 
galaxy is HI dominated (defui- -0.43). 

The radial profiles of dust surface mass densities are shown 
as green triangles in Fig. 
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The slope of the dust profiles tends 
to be intermediate between the molecular component (steeper 
slope) and the atomic component (flatter). 

To characterize the environmental effects in our sample, we 
show in Fig. 12 the bin-averaged profiles of the dust, atomic, 



molecular, and total gas surface mass densities for the three sub- 
samples defined above. In each panel red circles, black triangles, 
and blue stars, indicate the non-deficient, slightly Hi-deficient, 
and Hi-deficient subsamples. For the molecular component we 
show only the case for the Galactic Xco, since the conclusions 
about the environmental effects are similar even if the Xco de- 
fined in Eq.[T|is used. 

Galaxies with higher HI deficiency have steeper gas pro- 
files and lower surface mass densities, indicating that the lack 
of atomic hydrogen affects the formation of the molecular gas, 
and thus the total gas component. However this trend is not com- 
pletely obvious because, as shown by Fuma galli et al.| ( |2~009| l, HI 
deficiency is a condition necessary but not sufficient for molec- 
ular gas deficiency. In our subsample of Hi-deficient galaxies, 
according to the definition of H2-deficiency of Fumagalli et al. 
(2009), only NGC 4579 is strongly H 2 -deficient, and for this rea- 
son the trend observed for the atomic hydrogen in the bottom 
right panel of Fig. [12] is less evident for the molecular gas com- 
ponent (top left panel in Fig. 12 1. 



As the HI deficiency increases, the surface mass densities are 
lower for all the ISM components. However, for the atomic hy- 
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Fig. 9. Same as Fig.^for the five positions observed in NGC 4424. 



drogen in the central regions of the disks of Hi-deficient galax- 
ies, the surface mass density appears somewhat higher than that 
of slightly Hi-deficient galaxies. This effect could be understood 
as ram-pressure-induced compression of the gas in the inner re- 

200911 that increases 



gion of the disk (|Tonnesen & Bryan 



cally the gas density ( Byrd & Valtonen 



lo- 



1990). Alternatively, the 



differences in surface mass density can be simply due to uncer- 
tainties in the inclination corrections for two of the Hi-deficient 
galaxies, NGC 4388 and NGC 4402, which are nearly edge-on. 

Dust radial profiles of Hi-deficient galaxies (blue stars) de- 
crease steeply at R/R25 > 0.4, with respect to the non-deficient 
galaxies (red circles), but have comparable radial slopes in the 
inner parts of the disk, since the dust is more confined to the disk 
where ram pressure is less efficient in the ISM removal. Such a 
trend has already been found by Cortese et al. ( 2010| l. 



4.2. Dust-to-gas ratios 

In Fig. 13 we show the radial profiles of the dust-to-gas ratio 
for the three subsamples of galaxies, obtained after convolving 
the gas components to the resolution of the dust mass maps. The 



black solid line in each panel shows the bin-averaged profiles for 
each subset. 



In the top panels of Fig. 13 we use the constant Galactic 



Xqo- The averaged profiles show only moderate change with HI 
deficiency; the dust-to-gas ratio is almost constant throughout 
the disk for non-deficient galaxies (top left panel of Fig. 13 I, 



while it increases with radius for the strongly Hi-deficient subset 
(top right panel of Fig. 13 1. This could be explained because the 
gas, particularly the atomic component, is more easily stripped 
with respect to the dust ( [Cortese et al.| |2012[ ), thus producing a 
rapid increase of the dust-to-gas ratio at R/R25 > 0.2. However, 
the average slope of this subset may be biased by the highly HI- 
deficient galaxy NGC 4569 (defm= 1 .47), which has a gas profile 
that decreases rapidly at small radii. 

Similar trends with HI deficiency are found when the 
metallicity dependent Xqq is used (bottom panels of Fig. 13 I. 
However, there is a systematic change in the dust-to-gas ra- 
tio. Because of the lower Xqo in the metal-rich nuclei, the H2- 
dominated total gas mass surface densities decrease. Thus, the 
central dust-to-gas ratio increases, consequently producing neg- 
ative radial gradients. As shown in the bottom left panel of Fig. 
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Fig. 11. Dust (green triangles), atomic gas (red stars), and molecular gas (blue circles) surface mass densities as a function of R/R25 
for the extended sample. The molecular component has been estimated using the Galactic Xco (solid line) and the radially varying 
Xqo (Eq. T[ dashed line). The total gas mass surface density (for the Galactic Xco only), is indicated with black diamonds. The 
horizontal dotted line indicates £ = 9 M Q pc~ 2 . 



13 this effect is more evident for non-deficient objects, confirm 
ing what was already found by Magrini e t al.| ( |201 1) and by 
Foyle et ah] (|2012|l on the nearby galaxy M83. Also, Munoz- 



Mate os et ai.| ( 2009] >, using Spitzer data, find decreasing dust-to- 



gas profiles for an Xco factor that depends on the metallicity 
gradient. Simple considerations imply that the dust-to-gas ra- 
tio should follow linearly the metallicity; both dust-to-gas and 
metallicity trace metals, the former in the solid and the latter 
in the gaseous phase (Drain e~et al.||2007) . Thus, a lower-than- 
Galactic Xco in a galaxy's center appears to be necessary for 
the dust-to-gas ratio to have the same gradient observed for the 
metallicity. However, because of the large uncertainties in the 



absolute calibration of the metallicity indicators (e.g. Magrini et 
[ah] |2011| l, this statement is inconclusive. Indeed, independently 
of the metallicity, when the surface mass density is dominated 
by the atomic gas at all radii the dust-to-gas ratio decreases with 
radius, similarly to the metallicity gradient (see e.g. Bendo et al. 
|2010b| for NGC 2403, and |Fritz et al.|2012|[Smith et al.|2012| for 
the Andromeda galaxy). 



5. The linear relation between I, gas and F250: a pixel 
by pixel analysis 

Fig. [14] shows the surface mass density of H2 (top panel, using 
a Galactic Xco), HI (middle panel), and the total gas (bottom 
panel) as a function of the surface brightness at 250 fim images 
(F250) on a pixel-by-pixel basis. We chose this band because of 
its good sensitivity and better spatial resolution. Moreover, the 
surface brightness at one wavelength is preferred over the dust 
surface mass density because it is related to directly observable 
quantities; the scatter in the relation increases if the dust surface 
mass density is used instead of the 250 /mi surface brightness. 
Because of the poorer resolution of some of the HI maps, we 
smoothed all the data to a beam size of 36", and regridded all 
the maps to the same pixel size. Each pixel covers an area of 
~ 9 kpc 2 , except for NGC 4189 (at a distance of 32 Mpc), for 
which the area corresponds to ~ 30 kpc 2 . The top panel of Fig. 
14 shows that E# a is proportional to the surface brightness at 250 
fim. Considering all pixels above the 250 /mi surface brightness 
3cr threshold value of 0.07 Jy beam" 1 and 2 H: > 3cr = 1.15 M n 



2 , we found a best linear fit: 



log Z H2 = 1.23 ± 0.05 • log F 25 o + 1-59 ± 0.03, 
with a Pearson correlation coefficient r = 0.86. 



(2) 
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Fig. 13. Dust-to-gas ratio as a function of R/R25 for non-deficient (red circles), def m < 0.4, slightly Hi-deficient 0.4 < def m < 0.7, 
and Hi-deficient galaxies (blue stars), defm> 0.7. The solid line shows the profile for the new mapped galaxies. The black solid line 
shows the average profiles. The molecular content has been estimated using both the Galactic Xco (top row) and the Xqo of Eq.[TJ 
(bottom row). 

Table 4. Dust and gas properties in the sample 



Galaxy 




M m 


M H2 


M dlls ,/(M H1 +M H2 ) 






10 7 M G 


10 s M Q 


10 s M Q 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


NGC 4189 


4.47± 0.78 


22.7+0.7 


17.4+1.34 


0.011 


0.43 


NGC 4294 


1.15+0.23 


18.4+1 


0.73+0.2 


0.006 


0.04 


NGC 4298 


2.29+0.4 


5.6+0.9 


10.4+0.8 


0.014 


0.65 


NGC 4299 


0.52+0.1 


12.4+0.5 


1.1+0.2 


0.004 


0.08 


NGC 4351 


0.4±0.07 


3.4+0.4 


0.71+0.2 


0.01 


0.17 


NGC 4388 


1.2±0.24 


4.15+2 


7.13+0.6 


0.01 


0.63 


NGC 4424 


0.59+0.1 


3.97+0.2 


4.2+1.2 


0.007 


0.51 



Column 1: galaxy name; col. 2: dust mass (Da vies et al.| [2010 ); col. 3: integrated HI mass (Chung et al. 2009); col. 4: integrated H2 mass (this 
work); col. 5: dust-to-gas ratio; col. 6: molecular-to-total gas mass fraction. The molecular masses have been estimated using the Galactic X C o- 



The relation is tighter for non-deficient galaxies (red cir- 
cles) with respect to galaxies with a pronounced Hi-deficiency 
(blue stars). When computing a best linear fit we found a slope 
m — 1.29 + 0.07 (r = 0.91) for the non-deficient subsample, and 
m = 1.08 ± 0.1 (r = 0.81) for the Hi-deficient one. For both 
subsamples we find a similar intercept, ~ 1.6. The slightly bet- 



ter correlation between the 250 fim surface brightness and 
for non-deficient objects could indicate that the environment per- 
turbs both the molecular and dust disk, but in somewhat different 
ways, thus increasing the scatter. 



In the middle panel of Fig. 14 no clear correlation can be 
seen between Y,hi and F250. The spread is larger than for 
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0.5 

R/R25 



Fig. 12. Average radial surface mass densities of the dust (bot- 
tom left panel), Ehi (bottom right panel), (top left panel), and 
2 Hl01 (top right panel). Red circles, black triangles, and blue stars 
represent the average profiles for non-deficient (defm < 0.4), 
slightly Hi-deficient (0.4 < def m < 0.7), and Hi-deficient (def m > 
0.7) galaxies, respectively. The molecular component has been 
estimated using the Galactic Xqo- 



Mp, pc 2 for F250 ~ 0.3 Jy beam ') than Hi-deficient ones (~ 0.5 
"~ 2 ), though both samples contain outliers. This indicates 



and it appears that non-deficient galaxies have higher (- 

1 G pc~ 2 for F 25 o 
M G pc" 

that uncertainties in the defm parameter prevent a clear charac- 
terization of the environment effects on all galactic properties. 

Adding the HI contribution to all pixels with detected molec- 
ular gas, as shown in the bottom panel of Fig. 14 we find: 



log T, gas = 0.99 ± 0.04 • log F 250 + 1.63 ± 0.02. (3) 

The slope in this case is shallower than that in Eq.[2]and roughly 
linear. However, the goodness of the correlation does not change 
significantly when adding HI (r = 0.87). This is in contrast with 
the correlation found by Corbelli et al.| ( |2012[ ) between the in- 
tegrated properties: the total gas mass correlates better globally 
with the dust mass than the single gas components. A similar dis- 
crepancy is seen in the relations between star formation rate and 
gas mass surface densities; while a strong correlation is found 
between the global star formation rates and the total gas mass 
surface densities (Kennicutt, 1989), a resolved analysis gives a 
better correlation with the molecular gas surface density (Bigiel 
[eTaLi [20081 [Schruba et al.||20TT] l. 

If we use an Xqo that varies radially according to Eq. [T] the 
correlation between E gflJ and F250 is slightly degraded (r = 0.8), 
and the slope becomes shallower (m = 0.75+0.04). The radial 
variation of the Xqo produces in fact flatter radial profiles for the 
molecular surface mass density, which explains the shallower 
slope. 



monochromatic luminosity at 250 fim and the total mass of hy- 
drogen (using the Galactic Xco and calibrating the relation on a 
few nearby galaxies observed with Herschel and with gas mass 
estimates). With the dust temperatures measured on our maps 
(T ~ 20 K, see Da vies et al.||201 1) , their correlation can be con- 
verted to log T. gas = log F250 + 1 -62, in good agreement with Eq. 
[3] However, we had to assume that the 250 /urn luminosity arises 
from the same region occupied by gas, which is guaranteed for 
the molecular gas component only. 

6. Conclusions 

In this paper we present new IRAM-30m 12 CO(1-0) and 12 CO(2- 
1) observations of a sample of spiral galaxies in the Virgo cluster 
which are in the HeViCS field and in the VIVA sample. We ob- 
served four galaxies in OTF mode and three in PS mode. For ob- 
servations done in OTF mode CO-maps are discussed; for galax- 
ies observed in PS mode we only estimate the total molecular 
content. We include in o ur analysis a samp le of eight galaxies 
observed in 12 CO(1-0) by Kuno et al. ( |2007| , and investigate the 
effects of the environment on the surface mass density of the 
atomic and molecular gas phases and on the dust. Maps of the 
dust surface mass density have been estimated using the HeViCS 
data set (Davies et al. 2010 Auld et al. 2012[ ). Finally we an- 
alyze the relation between the different gas components and the 
dust distributions in relation to the defm parameter. Our main 
conclusions are the following: 

1 . For the new observed galaxies the global molecular gas-to- 
total mass fraction in the sample varies in the range 0.04 < 
fmoi ^ 0.65, using the Galactic Xqo, with the undisturbed 
galaxies showing the lowest values. Disturbed galaxies have 
a lower atomic and molecular content but the environment is 
more efficient at removing the HI than the molecular gas. 

2. The total gas surface mass density in the central regions 
of the disk, above a threshold of 9 M Q pc~ 2 is preferen- 
tially found in molecular form, in agreement with the results 
of Bi giel et"aT] ( |2008| l. However, a few Hi-deficient galax- 
ies show an order of magnitude lower threshold, indicating 
that the value may depend on local conditions. Highly HI- 
deficient galaxies have H2, HI, and dust surface mass densi- 
ties in the the central region lower by about a factor 2, sug- 
gesting gas and dust removal even in the inner disks. 

3. Galaxies with high HI deficiency have steeper radial profiles 
of the total gas, as a consequence of the combined effect 
of HI and H2 deficiency, in agreement with previous studies 
(Fumagalli et al. , 2009 ; Cortese et al. 2012[ l. The average ra- 
dial profiles of the dust surface mass density in Hi-deficient 
galaxies is steeper than in Hi-normal galaxies. This is con- 
sistent with the results found by Cortese et al. (2010 1, on the 
radial distribution of the flux density at 250, 350 and 500 /mi. 

4. The dust-to-gas ratio for non-deficient galaxies has a radial 
trend which depends on the Xco conversion factor, as already 
discussed by Magrini et al. (20lT). This ratio decreases ra- 



Eales et al. ( 20 1 2 1 proposed that dust emission can be used 



to map the ISM in high redshift galaxies: adopting the emissiv- 
ity of Galactic dust, they derived a linear correlation between the 



dially only if Xco increases radially due to metallicity gradi- 
ents. For Hi-deficient galaxies instead, the dust-to-gas ratio 
stays constant or increases radially, independently of the Xco 
factors used. This indicates that atomic gas is stripped more 
efficiently than the dust in a cluster environment. Since both 
dust and gas radial profiles truncate sharply well inside the 
optical disk, the dust-to-gas ratio cannot be traced at large 
galactocentric radius in highly perturbed galaxies. 
We observe a tight pixel-by-pixel correlation between the 
mass surface density of the molecular hydrogen and the 250 
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Fig. 14. Surface mass densities of molecular hydrogen (top 
panel), atomic hydrogen (middle panel), and total gas (bot- 
tom panel) as a function of the surface brightness at 250 fim 
on a pixel-by-pixel basis. Different symbols distinguish non- 
deficient (red circles), slightly Hi-deficient (black triangles), and 
Hi-deficient galaxies (blue stars). The pixel size has been set to 
36". The horizontal dotted lines show the 3<x level for the molec- 
ular and total gas (3<x =1.15 M Q pc~ 2 ). The diagonal lines show 
the best linear fits. Y, Hl in the top and bottom panels has been 
estimated using the constant Galactic Xco- A typical error bar is 
shown in the bottom panel. 



jj.m surface brightness. In Hi-deficient galaxies the correla- 
tion weakens, suggesting that the environment is perturbing 
both gas and dust. Adding the atomic gas to the molecu- 
lar component makes the correlation with the 250 //m sur- 
face brightness linear, but the correlation coefficient does not 
change significantly. Thus it is not clear if the dust emission 
at long FIR wavelengths can be used as a tracer of the total 
gas surface mass density, or of the molecular gas only. 
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Appendix A: Notes on individual objects 

A.1. NGC4189 

NGC 4189 is a relatively undisturbed Sc late-type galaxy (defui 
= 0.25) located at the edge of the Virgo cluster. The determina- 
tion of its distance and the membership to the cluster is debated, 
because Virgo is a cluster that is still forming and comprises dis- 
tinct clouds of galaxies at different distances. Using the relation 
1977]) in B-band, |Yasuda et al.| ([1997]) es- 



of Tully & Fisher 



timated for NGC 4189 a distance of 33.8 Mpc, i.e. more than 
twice the distance of the Milky Way from the center of the clus- 
ter (~ 17 Mpc). Solanes et al. ( 2002| l, starting from the HI dis- 
tribution, reconstructed the structure of the cluster, and found 
galaxies with elongated gas-tails, due to environmental effects, 
out to cluster centre distances of 25-30 Mpc. These results im- 
ply that at the distance of NGC 4189 from the cluster center it is 
still possible to find galaxies that are Virgo members. In the fol- 
lowing we considered for NGC 4189 the average distance of 32 
Mpc for the sub-group to which NGC 4189 belongs, determined 
using the Tully-Fisher relation in H-band ( Gavazzi et al. 1999 ). 
The HI disk is more extended than the optical disk, showing a 
symmetric warp (Chung et al. 2009) with enhanced emission 
in the South-East regions, in correspondence to a ridge of Ha 



emission (Koopmann & Kenney 20041. 



A.2. NGC 4298 

NGC 4298 is an Sc galaxy (def m = 0.41) located in the North- 
East region of Virgo at a projected distance of 3.2° from M87. 
The galaxy has flocculent spiral structure (Elmegreen et al. 
1999) and has a companion, NGC 4302, at a projected distance 
of 1 1 kpc, having a comparable velocity (Av « 30 km s~'). Both 
galaxies show Hi-tails pointing away from M87, probably due 
to ram pressure stripping and/or tidal interaction. However the 
undisturbed stellar component of NGC 4302 favors the ram pres- 
sure hypothesis (see |Chung et al. |2007| ). The HI disk in NGC 
4298 is compressed in the South-East region, in the direction of 
M87, as a consequence of the ram pressure exerted by the intra- 
cluster medium. The Ha emission is mostly concentrated in the 
central region of the disk, except for an enhanced arc in South- 
East direction, in correspondence with the HI compressed region 
( |Koopmann & Kenney) |2004[ ). 
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A.3. NGC 4299 & NGC 4294 

NGC 4299 is an Scd galaxy at a distance of 2.5° from M87 (defm 
= -0.43) which forms a possible interacting pair with NGC 4294 
{defm — -0.1 1). However both galaxies retain their HI gas which 
forms an extended HI disk ( Chun g~et al.| [20 09 ) with a long HI 
tail to the South West. Koopmann & Kenney (2004) showed that 
this pair has enhanced Ha emission in the northern part of the 
disk, where the HI is compressed. This indicates a stronger star 
formation rate, probably due to tidal interaction. 

A.4. NGC 4351 

NGC 435 1 is an Sc galaxy which is close to the cluster centre, 
about 1.7° from M87 and is slightly Hi-deficient (def m = 0.23). 
It shows a compression in the HI emission in the north eastern 
region, due to ram pressure ( |Chung et al.| |2009). The brightest 
optical region is offset with respect to the center, and also the Ha 
emission is asymmetric with enhanced emission in a circumnu- 
clear region (Koopmann & Kenney 2004| l. 



A.5. NGC 4388 

NGC 4388 is an Sab galaxy located at a projected distance of 
1.3° from M87. The galaxy is one of the most Hi-deficient in our 
sample (defm = 1.16), showing HI distribution truncated within 
the optical disk ( [Chung et alTl[2009] l. [Wshida et al.| ( [2002| found 
that the Ha emission extends for 35 kpc in the north east direc- 
tion. Vollmer & Huchtmeier ( |2003[ l found HI gas out to 20 kpc 
along the same direction as the Ha emission. Further observa- 
tions showed that the HI tail is even more extended, 110 kpc , 
with a mass of 3.4 x 10 8 M Q (jOosterloo & van Gorkom||2005). 



The most probable origin of this elongated tail is a ram pressure 
stripping episode that occurred about 200 Myr ago (Pappalardo 

|eTaTi[20T0l i. 

A.6. NGC 4424 



NGC 4424 is an Sa galaxy with def m = 0.97. |Chung etaLl(p007]l 
showed a long HI tail pointing away from M87 and Kenney et al. 
1996) observed star forming regions in the central kpc, but no 



star formation has been detected beyond these regions. Cortes 
|et al.| ( |2006| ) found two lobes in the molecular gas component 
in the central 3 kpc, with total mass M Hl ~ 4.0 x 10 7 M . The 
central regions are highly disturbed and asymmetric in the radio 
continuum, probably as a consequence of strong gravitational 
interaction ( Chung et al. 2009[ ). 
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